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1. INTRODUCTION 

Since the discovery of dark energy 0, 01 most of the attempts to understand it have consisted of guessing 
an equation of state as a function of the redshift parameter, z, evaluating one integral over z to determine the 
dark energy density and then another, necessarily numerical, integral to find a luminosity distance which can be 
compared with the experimental data Q . Thus it is not surprising that there is a lot of degeneracy - different 
equations of state give similar or even almost identical fits. What we do here is guess a form for the luminosity 
distance and confront our guess with the data straight away. Then we need only differentiate it once with respect 
to z to find the dark energy density and once more to find the equation of state. Thus we can, in principle, have 
analytic expressions for everything. 

This is not totally elementary, however, since a best fit to our guess for the luminosity distance, Dl{z), 
can result in an expression which violates some known property of Dl{z). These properties are automatically 
satisfied in the bottom-up approach. For example, the supernova data are bounded from above by the curve for 
dark energy alone, fix ~ ^tDl{z) ~ z{l + z), and from below by the curve for only matter, = 1,Dl{z) = 
2\J\ + z\\J\ + z — 1], so it natural to try a fit of the form 

d + azY 

U + az^ _ (1) 

an 

The best fit to this, using the data of Ref.|3| and a flat universe, is a = 4.61 and n ~ 0.534±0.004. Unfortunately, 
this grows faster than z at large z and, as we show below, that is not allowed because the density of dark energy 
would be negative for some values of z. 

The direct or top-down approach discussed here was introduced in a paper by Saini, Raychaudhury, Sahni and 
Starobinsky |^ for a flat universe. A reexamination of the results in this paper is included in what follows. In 
the next section we give some formalism and derive some conditions that the luminosity distance must satisfy. 
In Section 3 we make a few guesses for Dl{z), including the model of Ref. and flt them to the data. Finally 
in Sec. 4 we offer some discussion. 



2. FORMALISM 

The relevant parameters are the matter density and the curvature defined as fJjv/ = SttCpo/SHq and 
rife = —k/{R^H^) where po,Ho, and Rq are the current values of the matter density, the Hubble parame- 
ter, and the scale, and k = +1,0, or —1. With this method il^ is not a parameter which can be fit so, for 
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simplicity, we will fix Qm = 0.3. If flk is greater than or equal to zero then the (dimensionless) luminosity 
distance is given by 



Dr(z) 



1 + z 



sinh / dz 



where H{z) is the Hubble parameter 

H{z) 



= V^M (1 + + nx f{z) + f^fc (1 + Z)2 



with f[z) given in terms of the equation of state for dark energy, px = w{z) px, as 

dx 



f{z) = exp <^ 3 



1 + a; 



(1 + w{x)) 



and 



— 1 — Oaa — Ol, 



(2) 



(3) 



(4) 



(5) 



If flk is less than zero the the hyperbolic sine in Eq.© is replaced by the ordinary sine and f2fc is replaced by 
jrifcl = — rifc. However, after we take one derivative our expressions will be valid for all fi/j. Note that /(z) > 
for aU z and /(O) = 1. 

The simplest and most obvious of the conditions on D]^{z) is that it must be positive for all z. Also, for 
z << 1, D]^{z) must equal z, 



Dl(.z) > 0, 

Dl{z) = z, z « 1. 

Now to go further we need to solve Eq.® for /(z), so rewrite it as 



sinh ^ 



"ThDUz) 

1 + z 



/ rik Ho 



where the prime indicates a derivative with respect to z. Differentiating the sinh gives 



Dl{z) 
1 + z 



1 + z 



^il + zr + QkDliz) H{z) 



1 + z 



or, finally, using Eq. (Q, 

(1 + z)D'^{z) ~ Dl{z) ^ 
V(l + z)2 + rikDliz) ^QAiil + z)3 + (1 - Hm - ^k) f{z) + nk{l + zf 

Eq.lUni) is vahd for aU ilfc. 
Introducing D\{z) as 

/?i(z) = (l + z)i?;(z) - i?i(z). 



(6) 
(7) 



(8) 



(9) 



(10) 



(11) 



we note that, from Eq. (fTU]) . £'i(z) must be positive. This is not always easy to satisfy. Solving Ea. (fTn|) for /(z) 
gives 



(1 + Z)2 r(l + z)2 + 



D\{z) 



(1 + ^) - 



(12) 
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Since f{z) must be positive the terms in the { } brackets must be positive. On the other hand, for z » 1, we 
want the energy density in dark energy to be less than that of matter, i.e., 

(1 - riM - ^k) f{z) < nM (1 + zf. (13) 

So Eq.ini) gives 



nMii + z) + n,<^l±^l±^}^, auz, (14) 



and, from Ea. (|13|l we have 



(i + z)2 + nkDj(z) , , 

JJ^ [z) 

These conditions can be somewhat simphfied by writing Eq.Q for Dl{z) as 



.Tij: \ 1 + z 

where Dl{z) = Y{z) if flk is zero. Then f{z) of Eq. H12(l reduces to 



(i + z)2 r(i + z)2 



with 

Dy{z) = {1 + z)y'{z) -Y{z). (18) 

Di(z), given by Ea. ljlll) . is proportional to Dy{z) and thus Dy{z) must be positive. 
The constraint equations p4|l and (|15|l become 

0m(1 + 2) + Qk < ^-j^TV < 21]m(1 + ^) + rife, (19) 

where the first inequality holds for all z and the second is valid only for z >> 1. For z large the only solution 
which satisfies this and has Y{z) positive is 

Y{z) ^ Az ~ (20) 

with A and B positive constants. The A term, whose leading z behavior cancels in £'y(z), is necessary to keep 
Dy{z) positive. For the cosmological constant case, where f{z) = 1 for all z, B equals 2/y/TlM. For small z, 
Y{z) must, of course, go as z. Then, to leading order in z, Ea. ll9|l becomes 

nM < ^ < 2rjM (21) 

with 



Condition (|21|l is not perfect because it is only valid to leading order in y/z which is at most y/z = \/T089 = 33. 
In some models, such as Eq. (|29f) in the next section, the next term in Eq. H2(J|I . the constant term, is large and 
making it impossible to satisfy Eq. H19|) . Consequently, condition (|21|) is a necessary but not a sufhcient condition 
and for each fit to the data it is still necessary to check that the conditions in Eq. H19() arc satisfied. Eq. H21|l is 
only a guide as to where to look for a solution. 
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Once we find a model which satisfies these constraints, the dark energy equation of state is given by 

= -1+ + (23) 

which, using Eq. H17|) . can be written as 

^ 1 (1 + zrDyjz) + n.DUz) - 2(1 + zrViz) 
' 3Dy{z) {l + z)^-{l + z)nMD^z)-nkD^z) • ^ ^ 

Note that if we wish to satisfy the condition that w{z) > —1 then f'{z) > i.e., px must increase with z (just 
not as fast as pm). Wc wiU not worry about restricting w{z) in this paper. 



3. MODELS 



For the data we use the "gold" set of Ref. [j| plus the seven events given in Ref . 01 the total of which we will 
call the "essence" data. We also use the data set of Ref. which we call the "legacy" data. We compare these 
data sets with 5 logj^Q Dl{z) + M where M is an offset we fit as an extra parameter. WMAP |^ gives a comoving 
distance to the decoupling surface and wc translate this into a (dimcnsionless) luminosity distance and add it to 
whichever of the two data sets wc arc using Specifically, for this point we take Slogj^g Dl{z) = 17.79 ± 0.11 
at z = 1089 and we add it with no offset. 

From the previous section what is needed are expressions for Y{z) which give Dl{z) = z for z « 1, behave as 
Eq. H20() for z » 1, and have Dy{z) > for all z. The model which most closely resembles the form discussed 
in the Introduction but also satisfies these conditions is 



Y{z) = - Vl + az [VTTcz - 1] (25) 
c 



with 



C2 



e = - • (26) 
a 

For Qm = 0.3 the best fit to the essence data gives ^ = 0.6, a = 2.74, Qk = 0.112 and of 187.1 but, for 
large z, the ratio of energy densities px/pM is much too large, ~ 0.9. If we fix ^ at 0.35 a fit gives = 187.7, 
a = 2.43 ±0.20, flk = 0.029 ±0.036. In Fig.^wc show the 68.3% confidence contour in the a — ^ plane together 
with the best fit and the a — 2.43, ^ = 0.35 point. The latter gives a somewhat better ratio of energy densities 
at large z of 0.13, as shown in Fig.[3 For ^ below 0.35, px (/(z)) is not positive for all z. The legacy data with 
^ = 0.36 gives = 113.0, a = 2.50 ± 0.18 and ^Ik = 0.033 ± 0.035. The ratio of energy densities at large z 
is 0.16. Thus the two data sets give essentially the same answer for the parameters and results consistent with 
flk = 0. 

Another expression for Y{z) which satifies the condition (|20|) is 



Y{z) 

Here ^ of Eq. (|21|l is given by 



(27) 



4 /C\2n a^c^ 

— - 7 (28) 

Unlike the model given in Eq. I|25(l . the ratio of energy densities px/ Pm seems to be positive for all values of 
^ which satisfy condition 121|1 but, like the case of Eq. 125|) . this ratio is never very small. It turns out that the 
minimum occurs when a ~ c. Then to satisfy condition 121() together with Eq. H28() requires c < a, which, in 
turn, forces n to be very large in order to cancel the small denominator in Eq. (|28|) . This gives a lot of degeneracy. 
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Q. 1.0 




FIG. 1: The dashed curve is 68.3% confidence contour 
for the Y{z) of Eq. 12511 in the a — ^ using the essence 
data. Points in the shaded area satisfy the constraint 
< C ^ 2r2i\/. The cross is the minimum of a^nd 
the dot is the point a — 2.43, (, = 0.35. 



FIG. 2: The behavior of the ratio px{z) / pm{z) is 
shown for the Y{z) of Eq. 1251 using the essence data 
with ^ = 0.35 and Q.m ~ 0.3. 



Q. 1.0 




FIG. 3: The 68.3% confidence contour in the n — ^ 
plane is shown for the Y{z) of Eq. 1271 using the essence 
data. Points in the shaded area satisfy the constraint 
^ ^ ^ 2f2A/. The cross is the minimum of a^nd 
the dot is the point n — 0.5, ^ = 0.30. 



FIG. 4: The behavior of the ratio px{z)/pm{z) is 
shown for the Y{z) of Eq. 1271 using the essence data 
with 5 = 0.30 and ^ 0.3. 



For example, using the essence data, the minimum in occurs at 185.7 when the ^ parameter is 0.3, the other 
parameters are a = 0.918, c = 0.900, n = 34.57, flk = 0.104 and the ratio of energy densities at large z is 0.11. 
But if we allow to increase by one we can leave ^ unchanged, have n = ^ and a = 2.97, c = 0.351, fl^ = 0.141 
as shown in Fig.|3l In this case, the ratio of energy densities at high z is 0.17 and its behavior is illustrated in 
Fig-El The legacy data gives very similar results although allowing to change by one only gets n down to 0.9 
where a = 1.68 and c = 0.504. The only surprise is that the is larger than expected at 123.8. 

A brute force expression for Y{z) is 



Yiz) 



A' 



a + z 



2B 



b + z3/2 



(29) 
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with A, a, B, b, and c greater than zero and n greater than i. At large z the first and second terms on the 
right-hand-side reproduce Eq. H20() while at small z the last term gives the desired linear dependence on z. But 
this expression doesn't seem to work at all. A fit to the data gives seemingly reasonable values for and for the 
parameters; for example, using the essence data we find A = 4.80, a = 2.95, B = 1.8, b = 35.0, n = 1.92, c = 1.28 
and fife = —0.0883 with = 184.5. But the omitted constant term of Ea. (|20|l is equal to A{1 + a) which is large 
enough that the lower bound in Eq. H19|l is always violated for some values of z. 

Finally we consider the model given in Ref.0 which in our notation is 



Y{z) = 2(1 + z) 

with 



1 + z — ay^l + z — 1 + a 



+ z)+ -fV^ + z + 2-a-P-j, 



(30) 



e = ( ) ■ (31) 



For the essence data the minimum occurs at 



a = 1.93^°:?^ /3 = 0.66^°:°^ 

7 = 0-5lt"o"t = 0.09;°;J^ (32) 

with x^ — 184.8. This is smaller than those obtained for the models above, in particular model H25|l . but 
this probably isn't particularly significant; for fifc equal to zero x^ jumps to 188.4. The ratio of energy densities 
at large z is 0.12 which is about the same as model (|25|) . If we push the parameters to the lower bounds given 
by the errors the energy density ratio gets smaller but at the upper bounds it is unreasonably high ('--^ 1). The 
legacy data gives a reasonable fit with x^ = 112.5, but the central values of the parameters give a dark energy 
density which is negative for 2 < z < 40. Here, too, there are values of the parameters within the errors that 
result in a positive dark energy density for all z. 



4. DISCUSSION 



In Sec. 2 we derived some necessary conditions to be used as a guide in implementing this top-down approach 
for studying the dark energy. The simple models in Sec. 3 show that these conditions are not always sufficient. 
The model given by Eq. (|29|l satisfies the conditions but has no solution for the parameters which works for all z. 
The model given by Eq. I|27|) is worse in the sense that it has satisfactory solutions for an extremely wide range 
of values of the parameters. Knowing that a model does not work is useful since we can proceed to other forms. 
Having a model which works but where we cannot gainfully limit the parameters leaves us in a sort of limbo. In 
short the problem of degeneracy so prevalent in the bottom-up approach is no less severe here. In addition this 
model has a large bump in the ratio of the energy densities around z = 10 (see Fig.^J which seems unlikely to 
be physical. 

Thus we are left with the models of Ea. (|25|l or Ea. H30|) . For model (|25f) the best fit has an energy density 
ratio, px/pm, which does not go to zero at large z but wc can adjust the parameters to find a model where the 
ratio does get small without too much penalty in x^- This is the case shown in Figs.^andl^l Even then the 
model is not without problems. The equation of state w{z), as calculated from Eq. (|24|l . is less than —1 at z = 0, 
but we have agreed not to worry about conditions on w{z) here; even in a bottom-up approach it is not always 
possible to find fits to the data satisfying w{z) > —1 for all z. 

In some ways the model given by Eq. H30|) is the most satisfactory. It can have a small x^ if we allow fi^ to be 
nonzero and it can have small values for the ratio of dark energy to matter densities at large redshift if we allow 
the parameters to be one sigma off the minimum in the most favorable direction. However, if the parameters 
shift in the other direction the model rapidly becomes very unsatisfactory. Also here again w{0) is less than — 1, 
sometimes much less. 

On the positive side the data seem well behaved. For other reasons we think that Qk must be zero so it is 
pleasing to see that both data sets are consistent with this for the cases that wc tried. The predictions for the 



7 



parameters were also very similar for the two sets. However, a big part of this agreement may come from our 
inclusion of the WMAP point which has a strong normalizing effect. 

In summary it seems very hard to guess a form for the luminosity distance, Dl{z), which gives results that 
are physically sensible. This is not a problem in the bottom-up method. There it is all too easy to guess forms 
for w{z) which are sensible physically and automatically satisfy the conditions given above. Also the analytic 
solutions for the dark energy equation of state, which certainly exist in this top-down method, and are one of the 
main reasons for this approach, are likely to be rather complicated (see Eq. (|24|) '). So it is hard to know whether, 
with more data, this scheme of directly fitting the luminosity distances will be useful. 
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